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ABSTRACT 



The use of a self adaptive technique to compensate for the varying 
parameters of a flight control system is desirable for high performance 
aircraft and space vehicles. The technique discussed in this study con- 
sists of using a sinusoidal signal to measure the change in the parameters 
and of compensating for the change by a gain adjustment. Analog and digital 
computer methods are used to investigate the performance of the system and 
the accuracy of the Amplitude Closed Loop Response Criterion for measuring 
damping ratio. The effects of varying airframe flexibility, center of 
gravity, frequency of the sinusoidal signal, and the Reference signal value 
are included. The applicability of the self adaptive technique for a 
particular vehicle depends on t he dynamic characteristics of the airframe 
and on the desired na ximum variations of the performance characteristics. 

The writer wishes to express his appreciation for the assistance and 
encouragement given him by Dr. G. J« Thaler, Professor of Electrical 
Engineering at the U. S. Naval Postgraduate School. The writer also 
wishes to express appreciation to Mr. L. K. Mattingly who supervised the 
study during a field trip at Autonetics in the summer of 1961. 
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CHAPTER I 



INTRODUCTION 

The dynamic characteristics of an airframe vary when changes in 
speed, altitude, position of center or gravity, and dynamic pressure 
occur* The variation in the dynamic characteristics cause variations 
in the response of the airframe to commands from a pilot or program un- 
less some method of compensation is used. Many autopilots employ some 
form of gain variation for compensation. For propeller and early jet 
aircraft the gain variation was accomplished by a schedule which was a 
function of some air data parameter such as dynamic pressure, altitude, 
or Mach number. In space vehicles and high-performance supersonic air- 
craft the gain scheduling becomes very complex or impossible for the 
following reasons: 

1. The dynamic and static characteristics of the airframe change 
rapidly over wide ranges, 

2. Air data measurements are not accurate or are not available at 
high speeds and altitudes. 

3. The airframe dynamic characteristics are not knovm with suf- 
ficient certainty to permit scheduling. 

To eliminate the requirement for scheduling and air data measurements, 
automatic control systems are being developed which include a self adap- 
tive loop. The purpose of the self adaptive loop is to maintain optimum 
response of the airframe by compensating for the changes in the dynamic 
characteristics without the necessity for external measurements. 
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The Dither Self-Adaptive System is one method for obtaining optimum 
response by self compensation, A study of this system for pitch control 
was undertaken with the following objectives: 

1, To determine the ability of the system to correct for the 
varying dynamic characteristics of a flexible airframe, 

2, To determine the optimum variation of the damping ratio and 
natural frequency of the short period mode during a typical 
flight and how this variation is affected by airframe bending 
and by changes in the center of gravity of the airframe. 

The first objective was investigated by the comparison of digital and 
analog computer results. The second objective was Investigated primar- 
ily from an empirical approach although some analytical verification is 
included. 
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CHAPTER 2 



DITHER SELF ADAPTIVE CONTROL SYSTEM 
2, 1 Self Adaptive Systems 

Any closed loop control system has in effect some self-adaptive 
ability,, In the usual control system a desirable static output is obtain- 
ed under varying input and open loop parameter conditions by proper design 
of the feedback loops. The name **sel f-adapt ive” (I), however^ is reserved 
for only those control systems which maintain optimum dynamic response 
under varying conditions instead of the static response that normal feed- 
back optimalizes. 

Self-adaptive systems can be classified according to the ability to 
automatically compensate for either changes in the system input or changes 
in the system parameters^ such as environmental variations. The ideal self- 
adaptive system would compensate for both input and system parameter changes. 
Flight control systems are generally concerned with parameter changes and 
hence, are designed to compensate for only those changes. Since this in- 
vestigation is concerned with a flight control system^, only the ability of the 
system to optimalize the dynamic response under varying aerodynamic para- 
meters will be considered. 

There are three basic operations which an adaptive loop must perform: 

1) A continuous measurement of system dynamic performance - The 
dynamic performance measuring method should have a negligible 
effect on the system response signal and, in like manner, the 
system response signal should have negligible effect on the 
measuring method. 
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2) A coP-l inuous evaluation of the dynamic performance on the 
basis of some predetermined criterion - The selection of the 
evaluation criterion depends on the physical capabilities of 
the system. Factors^ such as realizability, complexity, and 
cost, may prevent incorporating into the system the ability to 
maintain a constant dynamic response. Compromises would have 
to be made so that the system would be physically capable of 
maintaining optimum a characteristic of the response; as, the 
relative stability - damping ratio - or the natural frequency. 

A figure of merit - which is a number or a method of compari- 
son - is then selected on the basis of accurately representing 
the selected dynamic response characteristic to be maintained 
constant . 

3) A continuous readjustment, based on the measured performance 
and evaluation, of system control parameters for optimum 
operation -- The methods used to adjust the parameters to 
obtain optimum operation can be listed under three headings: 
High Gain Linear Feedback, Programmed Compensation, and Compen- 
sation Using an Optimalizing Controller. 

In the High Gain Linear Feedback Method of adjusting performance under 
varying parameters the forward loop gain is maintained high by compensa- 
tion or with a relay. The desired dynamic resporse characteristics are 
contained in the feedback loop. The high gain with feedback makes the 
varying parameters in the forward loop negligible while the feedback loop 
determines the dynamic response of the system. This is not a new method and 
is used in many electronic applications for stabilization. 
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he Pr. gta-Tied ■ -i ^ t lor Method measures the varying parameter^ 
or the conditions whiv> vary the parameter^ and adjusts the compensating 
parameters in a controller in the forward loop to counteract the system 
variation. The adjustment is accomplished through a program which re- 
quires a knowledge of the relationship between the system and compensa- 
tion parameters. 

The Compensation Using an Optimalizing Controller Method is a type 
of feedback control which adjusts parameters so that a characteristic of 
the system response is optimalized; such as^ the damping ratio. In the 
Controller Method an error signal using the output response level of the 
s /stem can not be used as an error signal for the parameter adjustment as in 
a normal feedback method because the Controller is optimalizing a character- 
istic of the output and not necessarily the output level. The system may 
be used in conjunction with a model which gives the optimum characteristic 
in its output for a particular system input. The Controller samples the 
characteristic in the system output and compares the value with the optimum 
output of the model. The Controller then adjusts parameters until the model 
and system characteristic coincide, 

2.2 Description of Pitch Control System to be Investigated . 

The Dither Self-Adaptive System is capable of controlling a vehicle 
about its pitchy roll, and yaw axes. This investigation^ however, will 
be concerned only with pitch control. The system which is to be Investi- 
gated is shown in Figure 2.1. The position, rate and acceleration sensors 
that are required for the feedback loops are assumed to have a unity trans- 
fer function. The rectifier in the self adaptive loop is also assumed to 
have a unity transfer function. 
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In the convent i otial loop an error signal, S c, is formed by comparing 
the feedback signal to the command Input signal* The error signal is 
amplified by the adaptive gain, kg , and the resulting signal operates the 
servo actuator. This investigation involves a missile which has gimbaled 
thrust chambers. The servo actuator controls the direction of the chamber 
axis which results in a thrust vector signal <5Tt. The thrust vector signal 
affects the vehicle performance which is represented by the block labeled 
’‘Vehicle Dynamics”. The reaction of the vehicle airframe to the change in 
direction of the thrust axis is detected by position, rate, and acceleration 
sensors and is fed back through the Inverse Model to the input. The function 
of the Inverse Model will be explained in Section 2,3. 

For the self adaptive loop a dither signal is introduced at the input. 

This signal is then amplified by , operates the servo actuator, and 

« » 

excites the vehicle dynamics. The pitch acceleration, 0, is passed through 
a bandpass filter where the dither frequency is separated from the command 
input signal frequencies. The dither frequency is then rectified and com- 
pared with a reference signal. The ; operates an integrating 

servo which varies the adaptive gain, k ^ , until the reference and rectified 
signals are equal. 

2. 3 Principles Used to Maintain Optimum Dynamic Response . 

When a physical aerodynamic parameter changes, all the dynamic 
characteristic parameters of the system are affected. Typical physical 
aerodynamic parameters are Mach number and altitude, and typical dynamic 
characteristic parameters are natural frequencies and damping ratios. In 
terms of the system transfer function a change in aphy^'ical parameter could 
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shilc the posiilon of all the poles and 2 eros of the. transfer function* 

If constant dynamic response was desired for the system, the variation in 
all the poles and zero would have to be counter-balanced* The complexity 
of the situation can be visualized by placing a root locus plot alongside 
a missile traveling many times the speed of sound. The physical parameters 
of the missile are constantly changing, and each parameter would be varying 
all the poles and zeros on the root locus. A typical root locus with vary* 
ing poles and zeros is shown in Figure 2.2. A controller would have to be 
built to sense each of the variations and then decide how to counter* 
balance by compensation. In practice exact compensation is too complicated 
to be useful and approximations are always made. 

Simplicity is one of the main objectives of the dither adaptive loop. 

The only components that have been added to the existing ptich control 
system are a dither generator, filter, rectifier, and a comparison- integra* 
tion network. No variable shaping networks are used, and all adjustable 
compensation is accomplished by varying the gain of the system. 

Aviators are primarily affected physically by the short period mode 
of the pitch aerodynamic equations. Tests have indicated that pilots prefer 
a damping ratio of 0.7 and a natural frequency of 3 radians per second for 
this mode (2). Since the missile used in this investigation is designed 
for manned flights, the damping ratio of 0.7 and frequency of 3 were de* 
sired. Since the dither system does not use variable shaping networks, it is 
physically impossible to maintain both the damping ratio and frequency con- 
stant at the desired values. Therefore, a damping ratio of 0.7 with a 
minimum natural frequency variation are defined as the optimum dynamic 
characteristics of the short period mode during a flight. 
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There are actually two methods being employed in the overall system 
for varying the control parameter to obtain optimum performance* The 
control parameter of the system is the gain* The two methods^ which are 
discussed in Section 2.1, are High Gain Linear Feedback and Programmed 
Compensat ion. 

If the gain of the system could be maintained at a high value, the 

desired dynamic response characteristic could be inserted in a feedback 

loop and would determine the dynamic response of the overall system. 

This is the principle of the High Gain Linear Feedback Method shown in 

Figure 2.3. H(s) is referred to as the '^Inverse Model'* since the inverse 

of its transfer function determines the dynamic response of the system. 

If a damping ratio of 0,7 and a natural frequency of 3 radians per second 

2 

were desired, the Inverse Model would have a transfer function of s +4.2s+9* 
Unfortunately, the varying dynamic characteristics of the airframe, as 
represented by G2 (s), cannot be submerged by merely using high gain feed* 
back loops (3). The required aerodynamic control moments and power required 
for high gains in flight control systems are not available due to limits 
imposed by vehicle design considerations unrelated to the automatic control 
system design. Even if high gains could be obtained, the gains could re- 
sult in severe bending and slosh mode oscillations when considering flexible 
airframes. The values that the Inverse Model transfer function can assume 
and still affect the short period mode in the desired manner are also limited 
for flexible airframes. The limit on the values of the transfer function 
for this investigation is discussed in Section 2.5. 

For flight control systems the limitation on high gains and on values 
of the Inverse Model prevents the use of a true High Gain Linear Feedback 
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Pigure 2.5* General Feedback System. 
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Inverse Model can be used as a fixed 
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shaping nerwork to improve the dynamic characteristics of the short period 
mode and is employed as such in the system under invest igation* 

In view of restricting the Inverse Model to the idea of a fixed shaping 
network, the self adaptive scheme used in the Dither System can be thought 
of as a form of the Programmed Compensation Methods The three basic opera- 
tions that an adaptive loop must perform are listed in Section 2oi® The 
manner in which the Dither System performs these operations are listed below: 

1, The dither signal is the means for the continuous measurement 
of the system dynamic response^ In order to have a negligible 
effect on the system response signal, the dither frequency 
should be five or more times the short period mode frequency 
and the input level should be small so that the oscillations 
are not noticeable to the vehicle personnel* The reason for 
using pitch acceleration instead of pitch angle or rate as the 
pickoff for the adaptive loop is to enable smaller inputs for 
the same output* The upper limit for the dither frequency 
depends on the ability of the actuator to respond to the fre- 
quency and on the bandwidths of the rest of the system components. 
The upper limit is around 40 radians per second for available 
systems . 

2. The continuous evaluation of the dynamic performance is per- 
formed by comparing the amplitude of the pitch acceleration to 
dither input level evaluated at the dither frequency to the 
Performance Criterion Reference* It remains to be discussed 
below and in Chapter 4 how accurately the dynamic performance 
is being evaluated* In this investigation the damping ratio 
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is being used as the measure of dynamic performance* 

3. The continuous readjustment of the system control parameter, 
which is the gain, is accomplished by a servo driven by the 
error resulting from the evaluation procedure* 

In the Programmed Compensation Method the parameter adjustment is ac- 
complished through a program which requires a knowledge of the relationship 
between the system and compensation parameters* In the dither system the 
program consists of a single reference level which requires a knowledge 
of the variation of the system damping ration with the amplitude of the 
acceleration to dither input level ratio* 

The degree of self adaptability of the dither system depends on the 
accuracy with which the damping ratio is measured by the acceleration to 
input level ratio. There are many figures of merit that can be used as a 
criterion for the dynamic performance of a system. Some of the figures of 
merit are 

is an error signal. The applicability of these figures of merit for use 
in establishing a damping ratio of 0.7 has been determined for a second 
order system to be as follows: (4) 



selects zero damping ratio as optimum, 
selects 0.5 damping ratio as optimum. 

S selects 0.7 damping ratio as optimum, 

o 

selects zero damping ratio as optimum. 
selects 0.7 damping ratio as optimum. 

If the system was of second order, both ^ and 

would select 0.7 as the optimum damping ratio based on the figure of 

merit approaching a minimum value at a damping ratio of 0.7. The 

figure of merit has the sharpest minimum and would be the 
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preferred criterion. With a flexible airframe the system can not usually 
be approximated by a second order system with any accuracy 5 , and the ap- 
plicability of the above criterion for establishing optimum damping ratics 
of the short period mode under these conditions is not known. The corn- 



system is certainly more complex than the Dither System, and the possibility 
of the criterion selecting a damping ratio of 0.7 for a non second order 
system is considered doubtful. A different figure of merit called the 
Impluse Response Area Ration (IRAR) will permit the selection of any de- 
sired damping ratio and have the error signal assume a zero value when the 
response is at the selected damping ratio (5). The method required for use 
of the IRAR criterion would be to 'pulse the system and then measure the 
area of the response. The area is a measure of the damping ratio. By 
comparing the area with a reference area corresponding to the desired 
damping ratio an error signal is developed which would vary the gain of 
the system until the error was zero. In a flexible airframe control system 
which does not use shaping networks, the poles and zeros of the vehicle 
transfer function are continually varying. For the reference area to in- 
dicate the true damping ratio it would also have to vary with the parameters 
and in effect would need a self adaptive loop to adjust its value. The 
system would be very complex if not impossible to mechanize. 

The Dither System works on the principle of maintaining the gain of the 
system constant throughout the flight and assumes that a constant gain will 
maintain an optimum damping ratio. Referring to Figure 2.3 and substituting 



ponents required to mechanize the 




criterion into a 




At/Q an adaptive loop is formed. 
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has a denominator 

g 

of greater order than the rwimeratov ,-^{^.u/oJ= 

IVo is a large value. Therefore, is a direct measure of 

the gain, K^Ms j system. The self adaptive loop maintains the 

value of K;f equal to the reference value. If a constant value of 

is an accurate measure of the damping ratio^ ^ , for a flight 

control system, then a plot of HgMs versus ^ should show a constant 
^ particular for all flight cases. Using the I^S^S 

flight cases to be investigated a plot of KsMs vs cf is shown in 
Figure 2.4 for aerodynamic equations with no bending modes included, the 
first bending mode Included and the first and third bending modes included. 
The aerodynamic equations and conditions are discussed in Section 2.4. As 
is observed from the Figure, Ki is an accurate criterion for a con- 
stant for no bending, but as bending modes are included and the order 

of the system increases KgMf becomes a less accurate figure of merit. 

The investigation of the variation of the damping ratio for the various 
aerodynamic conditions and the determination of the optimum Performance 
Criterion Reference and the optimum dither frequency is contained in 
Chapter 4. 
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Aerodynamic Equations and Modifications . 

The airframe used in this investigation is the two staged Saturn 
Booster* Perturbated aerodynamic equations with coefficients for four 
flight cases are included in Appendix A* The flight cases are identified 
by the time in seconds after launch. The flight cases are t = 20^ t = 75, 
t = 101,9, and t = 200, The t = 200 case is for the second stage only, 
as the first stage separates before t = 200 seconds. The reference for 
the equations is the nose of the second stage. For flight cases t = 20, t = 
75, and t = 101,9 the acceleration, rate, and position sensors are located 
on the same platform in the first stage. The variation in performance for 
various platform positions will not be covered in this investigation. The 
platform position used had been selected for maximum system stability 
immediately after launch and for airframe considerations. The platform 
position for case t = 200 is in the nose of the second stage. The transfer 
functions which were determined from the equations using digital computers 
are included in Appendix A, 

The aerodynamic equations contain three bending modes and one fuel 
slosh mode, A typical root locus is shown in Figure 2,5, It will be 
observed from the Figure that the open loop system is unstable. Case t = 

20, which is the first case considered after launch and iv^here the missile 
is still verticle, is the only flight case with a stable open loop. An 
unstable aerodynamic condition occurs when the center of pressure is for- 
ward of the center of gravity, A missile of conventional shape will be 
unstable unless stabilizing fins are attached because the center of 
pressure of the fuselage alone is close to the shoulder whereas the center 
of gravity is close to the geometric center (6), For large missiles, such 
as the Saturn, the trend is towards stabilization through a control system 
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rather than providing for a stable airframe because of the weight required 
for the stabilizing fins. The control torques required for the control 
system are provided for in the Saturn by gimbaled thrust chambers. It 
will also be observed from the Figure that the Slosh Mode has a negligible 
effect on the system stability. Multitank design with slosh suppressors 
are used in the Saturn to obtain this condition. The different slosh mode 
frequencies correspond to different tanks. 

The second bending mode becomes unstable for an adaptive variable gain 
of approximately 0.1 for the t = 20, 75, and 101.9 cases. In order to 
eliminate the possibility of obscuring the actual dither system performance 
by the complexity of adding bending mode cancellation schemes, the unstable 
second bending mode was handled in one of two ways: 



compensation method was in use, such as a frequency tracker and notch 
filter scheme, which placed a pole and zero over the open loop zero and 
pole respectively of the bending mode. 

A bending cancellation method which could be used for procedure 1 



proposed by the Autonetics Division of North American Aviation, The purpose 
of the system is to form a signal of equal magnitude and opposite phase to 
the bending mode signal. The signal formed is then added to the control 
signal and effectively cancels the bending mode from the system response 
signal. Bandpass filters land 2 are variable tuned filters with center 
frequencies of Wo . If filter 1 is not tuned to the bending frequency. 



1. The second bending mode was eliminated from the aerodynamic 



equations assuming that a bending cancellation method was in use. 



2. The second bending mode was eliminated by assuming that a 



above is shown in ^ This is one of many cancellation methods 
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Figure 2.6. Bending Cancellation Method 



i:he bending frequen^^ coinpunent is phase shifted by an amount (p 

and the output signal from filter 1 Is sin (p^ o The bending sig-- 

nal is again phase shilled by an amount by filter 2 with output signal 
being sin^2^v2^-f The outputs from the two filters pass through 
phase shifters which produce a 90 degree phase shift between the signals. 
The shifted signals are then multiplied^ and the resulting signal drives 
an integrator which varies the center frequencies of the filters until the 
output from the multiplier is zero. As is shown below^ the output of the 
multiplier is zero when the filters are tuned to the bending frequency. 
Output of multiplier = sin(;^ 0 ^Z^^ cos {lUo. 1 ^ f Cp) 

= (sin cos (j) + cos sin^) (cos i^^.^coSf2^-sin I 4/0 t sin 2^). 

The integrator is not sensitive to frequencies much greater than the 
bending frequency. 

Therefore, 

Effective Output = 1/2 cos 2^ sin^ - 1/2 cos^ sin 2^ 

- 1/2 sin^ 

= 0 when ^ is zero. 

With the filters tuned to the bending frequency, the 180 degree phase shift 
is accomplished by subtraction from the feedback signal. The amplitude of 
the bending signal is adjusted by the rect i fier~ integrator comparison loop 
using the outputs of filters I and 2, for the comparison. Another bend- 
ing mode cancellation method using a rejection filter philosophy is con- 
tained in Reference 7. 

The different aerodynamic conditions that are investigated are listed 
below: 

1. No bending and no slosh modes included In the equations. 

The resulting equations are the rigid body equations. 
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2. Ihe iitsi bending mode added to the rigid body equationSo 

3. The first and third bending modes added to the rigid body 
equations* The second bending mode is assumed removed by a bending mode 
cancellation method* 

4. The first and third bending and first slosh modes added to 
the rigid body equations* The second bending mode is assumed removed by 
a frequency tracker and notch filter method* 

In Chapter 3 the analog computer investigation was conducted using 
aerodynamic conditions 1^ 2^ and 3 and flight cases t ~ 20^ 75^ 101*9 and 
200. The second bending mode was included in condition 3 for flight case 
t = 200 because the mode is not unstable for this case* 

In Chapter 4, the empirical investigation was conducted using aero- 
dynamic conditions 1,2,3, and 4 and flight cases t = 20, 75, and 101*9* 

2 . 5 Limitations on Inverse Model by Flexible Airframe * 

The system under investigation is unstable unless a shaping network 
is employed. From a root locus viewpoint if a complex pair of zeros were 
provided in the left half plane, the unstable short period mode would be 
shaped into the left half plane and the system would be stabilized* The 
Inverse Model provides these shaping zeros* 

If a damping ratio of 0,7 and a natural frequency of 3 radian per 
second were the desired dynamic characteristics of the short period mode, 
the Inverse Model should ideally provide a pair of zeros at s ^ -2*1 _ J2.18 
on the root locus. A high system gain would then establish a complex pair 
of roots of the system characteristic equation at values s - -2.1+j 2.18. 

The roots would fix the dynamic characteristics of the short period mode 
at a damping ratio of 0.7 and a frequency of 3* As was discussed in Section 
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2\3, bigh culpabilities eic 'iot available In flight eortrol systems. 

The desired dynar ic characteristics could still be. obtained even 
though high gains arc not available by Increasing the real and imaginary 
coordinates of the Inverse Model zeros. The zeros could shape the short 
period dynamic characteristic curve in such a manner that for lower values 
of gain the desired characteristics are available. Increasing the coordin- 
ates of the Inverse Model zeros is possible in a rigid airframe with 
negligible center of gravity movement. For flexible airframes the bend- 
ing modes must be considered in establishing the position of the Inverse 
Model zeros. 

If the coordinates of the zeros are Increased beyond certain values 
with a flexible airframe system, the zeros thape primarily a different 
mode than the short period mode. bliile the short period mode would be 
stabilized, its maximum damping ratio and natural frequency values would 
be less and its coupling with the shaped mode would be gi eater than when 
the zeros shape primarily the short period mode. A digital computer investi- 
gation of the optimum position of the zeros was conducted. 

The preliminary investigation showed that flight case t - 75 was the 
critical flight case and that the two bending modes plus slosh mode condi- 
tion was the critical condition. Flight case t - 75 was further Investi- 
gated as to the benefit of using rate and position feedback which would 
give one real axis zero or of using acceleration, rate^ and position feed- 
back which would give the complex pair of zeros. The complex zero arrange- 
ment provided for more shaping of the short period mode and less possible 
coupling between the first bending mode and short period m.ode. The coordin- 
ates of the complex zeros were then varied to obtain optimum values. 
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Typical root loci are shown in Figure 2o7. For Inverse Model position 3, 
the zeros shape primarily the first bending mode* While the short period 
mode is stabilized, it is not being shaped in an optimum manner and has 
the disadvantages of decreasing in damping ratio as gain increases beyond 
a certain value and of being coupled with the bending mode„ The short period 
mode could never obtain a 0,7 damping ratio. For Inverse Model positions 
1 and 2, the short period mode is being primarily shaped. The bending mode 
is also shaped further into the left half plane increasing its stability 
while producing a negligible increase in coupling. Position 2, with the 
zeros at s = 1.8 6, appears to be the optimum position. The maximum 

natural frequency of the short period mode at a damping ratio of 0,7 occurs 
at this position. As is observed, the natural frequency will never reach 
3 radians per second. A more complicated shaping network could be used to 
increase the natural frequency, but since one of the principal advantages 
of the Dither System is simplicity, a compromise has been made on the desired 
dynamic characteristics. As was defined in Section 2.3, the optimum dynamic 
characteristics of the short period mode for the investigation are a damp- 
ing ratio of 0.7 with a minimum natural frequency variation. 

The optimum Inverse Model transfer function for the severest flight 

2 

case and condition will remain fixed at s ^ 3, 6s+3.6 for all flight cases 
and for all aerodynamic conditions. 
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CHAPTLR 3 



COMPUIER SIMULATION 



3 ^ Int rodacl ion ^ 

The object of this Chapter is to investigate the ability of the 
Dither Self Adaptive System to correct for the varying dynamic character- 
istics of a flexible airframe during a typical flighto The correction is 
made by the adjustment of the control parameter which is the adaptive gain, 
kcT • 

An analog computer simulation was conducted for the primary purpose 
of determining the values of the adaptive gain which the system selected 
for the various flight cases and for the various aerodynamic conditions. 

The values of gain which were selected in the computer simulation are 
compared with the values of gain determined by digital computer methods. 
Digital computer techniques are employed in the investigation in Chapter 4, 
and the comparison with analog computer results in this Chapter is used as 
an indication of the accuracy of the results of Chapter 4 when the transi- 
tion is made from numerical to electrical simulation. The damping ratios 
of the short period modes which are associated with the gains selected by 
the analog simulation for the various aerodynamic conditions are determined 
from root loci using digital computer techniques and are compared to show 
the effect of a flexible airframe on the ability of the dither system to 
maintain a constant damping ratio. The damping ratios of the short period 
mode resulting from the use of the dither system are also compared with the 
damping ratios resulting from the use of a constant system. 

The dither frequency and reference level which would maintain through 
the flight an average damping ratio of the short period mode of 0.7 with a 
minimum variation of the damping ratio and natural frequency were not used 
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in the c..jputcr si-TL.'iC v»ri, icsults of the computer simulation 

indicated tbu r^ect^ssity rt determining optimum frequencies and reference 
levels. The investigation of this area is contained in Chapter 4. 

The secondary purpose of the analog computer simulation was to deter* 
mine for the various flight cases and aerodynamic conditions the rate at 
which the adaptive gain was adjusted and the reaction of the self adaptive 
loop CO wind gusts^ actuator movements^ and loss of the dither signal. A 
complete analysis of the seif adaptive loop was not conducted in the investi- 
gation; and therefore, the values of the loop gain and filter bandwidth for 
optimum dynamic response of the loop were not determined. The investigation 
was conducted using different loop gains for the various flight cases. A 
qualitative discussion of the adaptive loop is contained In Section 3«5. 

3. 2 Analog Computer Setup . 

The analog investigation was conducted using aerodynamic conditions 
1,2, and 3 as described in Section 2.4. For Flight Case t = 200 the 
second bending mode was included because the mode is not unstable for 
this case. The system was simulated in two, forty amplifier. Pace Computers 
by Electronic Association Inc. The voltmeter used to determine the K 
values and the integrating servo and rectifier used in the adaptive loop 
were inherent components of the computers. Two six channel Sanborn re- 
corders were used. A Hewlett-Packard, Low Frequency, Function Generator was 
used as the source of the dither signal. 

The computer simulation diagram is sho\^m in Figure B-1. The reference 
values of pitch acceleration, 0, pitch rate 0, and pitch angle, 0, which 

are measured at the nose of the second stage are converted to the platform 

« • » 

values, 9p , 0p , and Qp , for flight cases t ~ 20, 75, and 101.9, be- 
cause the sensors for these cases are on the same platform which is located 
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in the first stage. The reference values are the platform values for 
Flight Case t ~ 20 O 3 because the platform for this case is in the nose of 
the second stage. 

The first flight case simulated was t = 75 for aerodynamic condition 
3. Static and dynamic checks were conducted. The procedure used for the 
dynamic checks was to excite the bending and short period modes by voltage 
steps introduced into the amplifiers of the various mode simulations and 
to compare the resulting frequencies and damping ratios with the frequencies 
and damping ratios determined by digital computer techniques. The trans- 
fer function resulting from the removal of the poles and zeros associated 
with the second bending mode from the vehicle dynamic transfer function, 
which was derived from the aerodynamic equations with three bending modes 
and one slosh mode, was considered a good approximation for determining the 
first and third bending mode frequencies for condition 3. A comparison with 
the analog results of the dynamic checks indicated that new transfer func- 
tions based on the applicable aerodynamic equations were required. The 
transfer function for flight case t = 75 for aerodynamic condition 3 in- 
dicated that the only poles and zeros of the transfer function^ which was 
derived from the aerodynamic equations with three bending modes and one 
slosh mode, that were appreciably moved were the zeros associated with the 
third bending mode. The value of the zeros were changed from s =-“11 . 56+ j21. 
63 to s = “.17+; J24.3. The effect of the bending and slosh mode parameters 
on the vehicle dynamic characteristics can be determined by observing the 
movement of the poles and zeros of the transfer functions, which are contain 
ed in Appendix A, for the various aerodynamic conditions. The static checks 
agreed to within one half of a per cent with analytical values. The mode 
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frequencies of the analog computer agreed to within two per cent of the 
frequencies determined by digital methods^ and the damping ratios agreed 
to within five per cent. 

The procedure used to simulate the self adaptive loop was as follows: 

1. The filter was simulated and the bandwidth was varied until 
a sharp peaked frequency response was obtained with the peak at the dither 
frequency of 30 radians per second, A bandwidth of 3 radians per second 
produced a sharp peak, 

2. Two diodes which were inherent components of the computer 
were used to form a rectifier. The outputs from the two diodes were of 
different values; and therefore, gain pots were included with each of the 
diodes so that equal outputs could be obtained. 

3. The adaptive gain, kj , was simulated by using a linear 
gain pot driven by an integrating servo/ The servo and pot were inherent 
components of the computer. One volt into the servo resulted in a gain 
change of 0,2, A linear pot was used because it was readily available, 

A shaped gain pot should be used in reality to obtain a small gain to servo 
input ratio for low values of gain and a large ratio for large values of 
gain, A small gain to servo input ratio is desired at low values of gain 
so that a wind gust or other disturbance would not reduce the gain apprec - 
ably with the possibility of the system going unstable. As the gain in- 
creases the dynamics of the short period mode improve so that while a 
disturbance at high gains would vary the gain more than at low gains the 
response time is Improved and stability considerations are not as critical. 
Compromises must be made in the shaping of the gain pot between the large 
ratio of gain to servo input desired for fast adaptive response, the small 
ratio desired when disturbances are experienced, and the ratio required for 
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constant dynamic characteristics of the adaptive loop Itself during 
a flight. The discussion of maintaining the adaptive loop gain appreciably 
constant during a flight is contained in Section 3<,5* 

4. A gain pot was inserted before the integrating servo to 
provide for varying the ratio of to servo input and the adaptive loop 

gain. The adaptive loop gain was adjusted by varying the gain pot^ , 

until the gain adjustment response^ , ceased to oscillate and was 

critically damped vjith a suitable rate of change of the value of , 

The Reference DC signal was then adjusted to provide for a value of 
of 1. The adjustment of the adaptive loop gain to give a suitable rate of 
change of and the setting of the Reference level so that had a value 
of 1 was performed for each aerodynamic condition of Flight Case t ~ 75 
in order to obtain reference values. The reference values that were obtain- 
'ed eliminated the necessity for a complete evaluation of the analog simulat- 
ed loop and were used in the other flight cases to form the basis for the 
comparison of the values of that the system selected. 

When the system was energized, the servo chattered because it was 
following the half cycle sinusoidal like waveform at its input, A smooth- 
ing network would be required at the rectifier output to produce a flatter 
output and eliminate the servo chatter. In order to maintain the adaptive 
loop as simple as possible and also to investigate the effect of bending 
interference the bandwidth of the filter was increased to 5 radians. The 
third bending mode frequency which was now at the limits of the bandwidth 
was passing through the loop because the waveform from the rectifier flat- 
tened out so that the servo stopped chattering. Since a dither frequency 
of 30 radians per second excited the third bending mode for aerodynamic 
condition 3, the computer results would indicate the performance of the 



30 



( 

• I 



di t ^ 



lev?, when bending or noise signals were contained in the filter- 
ed fi cy of the self adaptive loop. In order to determine the per- 

formarc v r the system with no bending and minimum noise signals in the 
self adaptive loop, reference values x^ere also recorded for Flight Case 
t = 75 for aerodynamic condition 3 with a dither frequency of 20 radians 
per second. The filter was adjusted for a center frequency of 20 and a 
bandwidth of 3. The lox>;er frequency did not cause the servo to chatter. 

The general procedure used in the simulation for Cases t - 20, 101.9, 
and 200 ^^as as follox>7s: 

1, The computer was set up for condition 3 of the particular 
flight case. 

2, Static and dynamic checks were conducted. 

3. The Reference level and KsA x^ere set at the values used 
for the corresponding aerodynamic condition and dither frequency of Flight 
Case t -- 75. 

4, The values of as read from the computer voltmeter were 
recorded for aerodynamic conditions 1, 2, and 3 with a dither frequency 
of 30 radians per second and for condition 3 x^ith a dither frequency of 
20 radians per second. 

5. The value of KsA was adjusted to give a suitable rate of 

change of for the particular aerodynamic condition and flight case. 

The Reference level was then adjusted to give the corresponding value of 

as determined in Step 4 with a dither frequency of 30. 

6. Recordings were made at the dither frequency of 30 for 
aerodynamic conditions 1 and 3. During the run for each condition, 

was offset by a step voltage put into the integrating servo amplifier; an 
actuator signal, , \<fas put into the system by a step voltage into the 



31 







I 

I 



simulation for the actuator; a v;lnd gust was simulated by a step voltage 
into the attack angle, o( , amplifier; and a dither signal failure was 
simulated by turning off the signal generator. 

Step 6 was also conducted for Case t = 75. 

3. 3 Comparison of System Selected Gains . 

The values of the adaptive gain, , which were recorded by the 
procedures described in Section 3.2 during the analog computer simulation 
are compared in this Section with values of gain determined from digital 
methods to establish a correspondnece between the two procedures. A Recomp 
II Computer v;as used for the digital calculations. The damping ratios of 
the short period modes corresponding to the recorded adaptive gains are 
compared to establish the relative adapt ibility of the system under vary- 
ing aerodynamic conditions and dither frequencies and are compared with 
the damping ratios corresponding to a fixed gain system to establish the 
benefit of the adaptive system for the airframe being used for the investi- 
gation, Tire analog computer simulation was conducted for aerodynamic 
conditions 1, 2, and 3 as described in Section 2.4. The digital investiga- 
tion included aerodynamic condition 4, and the damping ratios associated 
with the gains for a particular Dither Reference level for this aerodynamic 
condition are included in the comparison of the damping ratios of the ana- 
log computer. 

The determination by digital methods of the values of which the 
system would select was accomplished by computing the amplitudes of the 
closed loop response for various gains and for frequencies of 20 and 30 
radians per second. Plots of the amplitudes of the closed loop response 
versus for a dither frequency of 30 were made of the flight cases for 
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each aerodynamic condition. For aerodynamic condition 3 plots were made 
for dither frequencies of 30 and 20, The plots are shown on Figures B-2 
(a) through B-2(f). An adaptive gain of 1 was set on the analog computer 
for each aerodynamic condition of Flight Case t = 75 and was used as the 
standard for determining the values of Reference level and adaptive loop 
gain pot settings which were maintained constant in recording the adaptive 
gains of the remaining flight cases. The adaptive gain value of 1 for Case 
t = 75 was also used as the standard for the digital method. The reference 
amplitude of the system was determined on the plots by the intersection of 
the equal to 1 line with the Flight Case t = 75 curve. The reference 
amplitude eliminated the necessity of evaluating the components of the 
adaptive loop to determine a Reference Signal level and was used to deter- 
mine the values of gain which the system would select for the other flight 
cases by the intersection of the reference amplitude and the particular 
flight case curve. 

The comparison of the adaptive gains determined by analog and digital 
methods is shown in Table I and II for aerodynamic conditions land 2 
respectively at a dither frequency of 30, The reference amplitude deter- 
mined for aerodynamic condition 3 at a dither frequency of 30 was too large 
for use with Flight Cases t = 101,9 & 200 as shown on Figures B-2(c), The 
dither frequency was varied in the digital method for each flight case 
until a reference amplitude was determined i^hich resulted in the gains which 
were selected by the analog computer. The dither frequencies required to 
obtain the gains are shoTO in Table III, The plots of closed loop response 
amplitude versus for the different dither frequencies required are shown 
in Figure B-2(d). The third bending mode was excited at a dither frequency 
of 30 radians per second and the filter bandwidth allowed the bending signal 
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TABLE I 



Gain Comparison for Aerodynamic Condition 1 at Dither Frequency of 30 



Flight Case 


Gain Selected by 
Analog Computer 


Gain Determined by 
Digital Procedure 


t = 20 


1,73 

1 


1.69 


t - 75 


1.0 


1.0 


t = 101.9 


0.54 


0.535 


t = 200 

( 


1.15 


1.17 



TABLE II 



Gain Comparison for Aerodynamic Condition 2 at Dither Frequency of 30 



Flight Case 


Gain Selected by 
Analog Computer 


Gain Determined by 
Digital Procedure 


rt 

II 

O 


1.35 


1.3 


t = 75 


1.0 


1.0 


t = 101.9 


; 0.8 


0.72 


t = 200 


3.4 


3.42 



3A 



to pass thrcagh the adaptive loop as was discussed in Section 3.2. 

The dither signal source could not be accurately set to within 0.628 
radians per second, and since two dither frequencies were used for each 
flight case, the signal source had to be adjusted for each case. The 
variation in the digital frequencies required to obtain the analog gains 
is approximately within the accuracy of the signal source settings. The 
necessity for very accurate determination and setting of the dither fre- 
quency and Reference Signal level when bendipg modes are excited close to 
the dither frequency is illustrated by the difficulty in obtaining a compari- 
son of the gains for the digital and analog methods. When bending or noise 
signals pass through the adaptive loop, there is a possibility of the system 
having a choice of selecting either of two adaptive gains which are relative- 
ly close in m.agnitude. The selection of the wrong gain might cause the 
system to be unstable. As is shown in Figure B-2(d), the reference ampli- 
tude intersects the flight cases at two relatively close values of 
The comparison of the gains determined for aerodynamic condition 3 at a 
dither frequency of 20 is shown in Table IV. The difficulties experienced 
with a dither frequency of 30 were not experienced at a frequency of 20 
where bending modes are not excited. The plots of closed loop response 
versus for a dither frequency of 20 are shown in Figure B-2(e). 

Root loci x«jere determined for the various flight cases and aerodynamic 
conditions using a Recornp II digital computer. The damping ratios of the 
short period modes corresponding to the gains selected by the analog computer 
with a dither frequency of 30 and to the gains selected by digital techniques 
for aerodynamic condition 4 v;ith a dither frequency of 30 were determined 
from the root loci and are compared in Table V. 
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TABLE III 



Frequencies Required by Digital Method for G^ins Determined by Analog 
Method for Aerodynamic Condition 3 at Dither Frequency of 30 



1 Flight Case 



Gain Selected by 
Analog Computer 


Frequency Required 
by Digital Method 
in Radians per Second 


o 

CM 

II 

U 


1.0 


28.5 


j t . 75 


1.0 


28.5 


t 

i t = 101.9 

i 


2. 1 


30.0 


\ j 

; t = 200 ; 


; 5.0 


ro 

CO 



TABLE IV 

Gain Comparison for Aerodynamic Condition 3 at Dither Frequency of 20 



Flight Case 


Gain Selected by Analog 
Computer 


Gain Determined by 
Digital Procedure 


= 20 


1.35 


1,42 


t = 75 


1,0 


1.0 


t = 101.9 

1 


0.6 

i 


0.6 


t = 200 ' 


1 

1.5 


1.5 

1 
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TABLE V 



Damping Ratio Comparison for Flights Under Various Aerodynamic Conditions 



Flight Case 


Damping Ratios 




Damping Ratios 


Damping Ratios 


Damping Ratios 




for 




for 


for 


for 




Aerodynamic 




Aerodynamic 


Ae rodynamic 


Aerodynamic 




Condition 1 




Condition 2 


Condition 3 


Condition 4 


t = 20 


.82 




.78 


.725 


.585 


t - 75 


.794 




.838 


.843 


.72 


t = 101.9 i 

1 


i .838 

1 




.87 


.912 


.875 


t = 200 ! 


1 

1 .89 

. _ _ _ . _ . 




,924 


.93 


.89 


Total 1 


i 


f 








Variation j 


i 








in Damping ! 


( 1 \ 


i 


.14 


,21 


.21 


R'' t io 




L 








Average 


1 








Damping 


\ 






,792 


Ratio for 


.836 


1 


.853 


.852 


Flight 




4 

t 

i 












t 








Variat ion 




I 








1 from t = 75 


i + . 09 : 


A 


+ .086 


+ .087 


i +.17 


1 

Reference ; 


I -.0 


9 

i 

L 


-.058 


-.118 


i -.035 


' 


i 
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Flight oa.st t 7'i js jsed as the reference for both the analog and 
digital procedures; and there tore^ the ir.axiirmm variations of the damping 
ratios of a particular flight from the damping ratio of cases t = 75 is 
probably the best criterion to show the effect of bending modes on the 
adaptibility of the system. The results of the comparison of the fipst 
three aerodynamic conditions shows adaptibility decreasing as bending is 
included, but the comparison is valid only for a dither frequency of 30 
and must be checked for other dither frequencies. The investigation as 
to the effect on the results of a comparison when the dither frequency is 
varied is contained in Chapter 4. A different method was used in aero- 
dynamic condition 4 than in condition 3 to remove the second bending mode. 

The different methods resulted in a large variation of the position of 
the zeros associated with the third bending mode as is noted by examining 
the transfer functions in Appendix A. The comparison of the two conditions 
in Table V indicates that for a dither frequency of 30 the frequency track- 
ing and notch filter scheme for removing the second bending mode results in 
a more adaptive system than the bending cancellation scheme. The total varia 
tion of the damping ratios of the two conditions is the same but condition 4 
includes a slosh mode, has the larger value of variation from the reference 
value of the t = 75 case towards greater damping, and has an average damping 
ratio during the flight less than condition 3. For a larger average damping 
ratio during a flight, a smaller variation in the damping ratios would be 
expected because a larger average damping ratio means larger gains are used 
and large gains force the damping ratios towards the constant damping ratio 
determined by the position of the Inverse Model. 

The damping ratios as determined from the root loci for the gains select 
ed by the analog computer for aerodynamic condition 3 with dither frequencies 
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of 30 ar-.i • 1.i lablt vi with the damping ratios of a fixed 

gain system v’ith the ii sot at a value of 1. The gain of a value of 1 
was selected for the fixed gain system in order to establish the same 
reference for Case t - 75 that was used for the Dither System. The results 
of the comparison show that the fixed gain system is more adaptive than 
the dither system at a frequency of 30 radians per second but is less 
adaptive than the dither system at a frequency of 20, The great difference 
in the performance of the dither system as the dither frequency was varied 
led to the investigation of obtaining the optimum dither frequencies for 
the different aerodynamic conditions. This investigation is contained in 
Chapter 4. 

3.4 Reaction of the System to Disturbances . 

Aerodynamic conditions 1 and 3 were used in the investigation of the 
rate at which the adaptive gain was adjusted and the reaction of the adap- 
tive loop to wind gusts, actuator movements, and loss of dither signal. 

The investigation was conducted with a dither frequency of 30 radians per 
second so that the evaulation of the adaptive loop for aerodynamic condition 
3 would be performed under the unfavorable conditions of having bending 
signals included with the filtered dither signal in the adaptive loop. 

Since a complete analysis of the adaptive loop was not performed, the 
adaptive loop gain was adjusted for each flight case with the following 
considerations : 

1. To produce a suitable rate of change of the adaptive gain, 

2. To show the variation in Che loop performance for different 
flight cases with the same setting of the servo gain pot, 

3. To use several values of for the flight cases with one 

value producing an overshoot in the response. 
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TABLE VI 



Darpii'g Ratio Coirparison of Dither System and Fixed Gain System 



Flight 1 

Case 1 


Damping Ratios 
for 

Dither Frequency 
of 30 


Damping Ratios 
for 

Dither Frequency 
of 20 


Damping Ratios 
for 

Fixed Gain System 


t = 20 


.725 


.78 


.725 


t = 75 


.843 


.838 


.843 


t = 101.9 


.912 


.844 


.88 


t = 200 


.93 


.9 


.882 


Total 
Variation 
in Damping 
Ratio 


.21 


. 12 


. 16 


Average 
Damping 
Ratio for 
Flight 


.852 


.84 


.832 


Variat ion 
from 
t = 75 
Reference 


1 

i +.09 

j -.118 


4.06 

-.06 


+ .039 
-.118 
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The ^ 1 £ V r corresponding flight cases under the 

difitjCc^ - ‘ i I : ■ It tons to show the effect of bending on the 

loop per c. 

Recordi':L:> were ujade of the reaction of the system to various distur- 
bances for the different flight cases and aerodynamic conditions. The traces 
for the flight cases for aerodynamic condition 3 showing the actuator 
command signal the attack angle^ 0( ; the adaptive gain response, ; 
the signal from the adaptive loop rectifier, SAE; and the first bending 
mode, dj, are contained in Figures B-3(a) through B-3(d), The procedure 
used in obtaining the traces is contained in Steps 5 and 6 of the procedure 

used to simulate Cases t - 20, 101,9, and 200 in Section 3.2, 

The comparison of the rates of change of when Kg is offset above 
and below the undisturbed magnitudes selected by the system for the vari- 
ous cases and conditions is shovjn in Table VII, Cases t - 20 and 75 have 
the same setting, or Kg to servo input ratio, but th.. response tim.es 
vary considerably. The variation illustrates the fact that the Kg to 
servo input ratio must also be programmed for the flight if constant dynamic 
conditions are to exist in the system. There must be an optimum Reference 
Level, Dither Frequency, and Adaptive Loop Gain for optimum operation of the 
dither system. The program required to maintain the optimum value of the 
Kg to servo input ratio can be placed into the loop in the form of a shaped 
gain curve for pot as is explained in Section 3.5, The adaptive loop 

gain for Case t 101.9 resulted in an overshoot in the Kg response as is 
shown on Figure B-3(c) The response time for this case is the fastest as 

would be expected since the Kg responses for the other cases are ov^i. damped. 

As bending was added the response time increased except for Case t - 101,9 
where the tesponse time remained constant. In order to form some basis 
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TABLE VII 



Coinpai : son -U Kattj uf Change of the Adaptive Gain for Various Loop Gains 

and Aerodynamic Conditions 



Flight 

Case 


Lindisturbed 
slalixes of 

H 


Time (Seconds) 
per Unity Change 
in 


Gain Setting in 
Adaptive Loop 


Aerodynamic 

Condition 




1 3 


1 3 


offset abc^’''] 
t = ^0 1 

offset belowj 
j 


\ -.0 

[ 

i 

1 


3.4 22.0 

3.8 4.5 


.025 .025 ; 

.025 .025 

i 


offset above 

t - 75 

offset below 


1.0 


5.0 14.0 

5.0 15.0 


.025 .025 ] 

.025 .025 { 


offset above 
t - 101,9 

offset below 


2.1 


2.6 2.6 

4.0 4.0 


0.2 0.2 

0.2 0.2 


offset above 

t = 200 

offset belo\ 


5.0 


6.6 11.0 

4.4 7.2 


0.0 0.1 1 

1 

0.0 0.1 



I 
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for establishing if the ti^;.as listed in Table VII are sufficiently fast for 
the system to function properly^ the time intervals between flight cases 
are divided by the gain changes required between flight cases to estab- 
lish figures of merit. The figures of merit are based on the assumption 
that is changing linearly with the time between flight cases and are 

listed below: 

t = 20 to 75 o^seconds per unitychange in 

t = 75 to 101.9 24,6 seconds per unity change in 

t = 101.9 to 200 33.8 seconds per unity change in 

Even though the loop was not optimized the times listed in Table VII are 
well within the values established as figures of merit. 

The comparison of the reaction of the adaptive loop to actuator signals 
for the various flight cases and aerodynamic conditions is shown in Table 
VIII. The magnitude of was changed in all cases and conditions when an 
actuator signal was introduced into the system^ but the adaptive loop re- 
turned , while the signal was still being applied, to the original 

magnitudes for all cases and conditions in relatively small time intervals. 
The ideal system would result in no change in when an actuator signal 
was experienced. Excessive per cent changes in occurred for Flight 
Case t = 20 for aerodynamic conditions 1 and 3 and for Flight Case t 
101.9 for aerodynamic condition 3. If the adaptive loop had an optimized 
gain with a properly shaped g^in pot instead of the linear pot used in 
the simulation, these large changes in v;ould be greatly reduced. One 
of the considerations that must be used in arriving at the proper shape for 
the gain curve of the pot is to have a small ratio of to servo input for 
small values of , The undisturbed magnitudes of Kg for Case t ^ 20 was 
1 and for t = 101.9 was 2.1 so that the shaped pot would aid these cases 
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Table viii 



Cotnpar i soi. ui the Reaction of the Adaptive Loop to 
Actuator Signals 



Flight Case 
Aerodynamic 
Condition 
3 


(Sc Step 
Input in 
Degrees 


Change in 
Angle of 
Attack^ 
in Degrees 


Change In 
K^- 


Per Cent 
Change in 
Ki- 

■ 


Time(Seconds) 
to return to 
Undisturbed 
Value of 


Ga In 
Setting 
in Adaptive 
Loop 


t = 20 


-0.6 


1.0 


-0.2 


20.0 


4.6 


.025 




+0.6 


1.0 


-0.2 


20.0 


4.8 


.025 


t = 75 


-0.4 


0.6 


-0.08 


8.0 


5.0 


.025 


1 


+0.4 


0.6 


+0. 16 


16.0 


20.1 


.025 


t = 101.9 ' 


-0.15 


0.6 


+0.8 


33.0 


4.8 


0.2 




+0.15 


0.6 


-0.8 


33.0 


3.0 


0.2 


t = 200 


-0.62 


0.5 


-0.4 


7.1 


9.0 


0.1 




+0.62 


0.5 


+0.4 


7.1 


8.5 


O.T 


. Flight Case 














Aerodynamic 














Condition 1 














t = 20 


-0.85 


1.1 


-0.2 i 


20.0 


3.0 


.025 




+0.85 


1.1 


+0.2 


20.0 


3.0 


.025 


t = 75 1 


-0.7 


0.7 


-0.11 


10.0 


3.6 


.025 




+0.7 i 


0.7 


+0.11 


10.0 


3.8 1 


.025 


t = 101.9 


-3.2 


1.0 


+0.16 


7.2 


3.0 


0.2 




+0.2 


1.0 


-0.08 


3.6 


2.0 


0.2 


t = 200 


-0.5 


0.45 


-0.2 


3.8 


8.0 


0.1 


1 1 


+0.5 

! 


0.45 


+0.2 


3.8 


7.8 


0.1 



in prevent LKg latgo in Kg when disturbances were experienced. 

The same magnitude actuator signals for the bending aerodynamic condition 
produced approximately the same change in the angle of attack as for the 
no bending aerodynamic condition^ but the bending condition resulted in 
greater changes in and longer time intervals to return to conditions 

that existed before the disturbance than the no bending condition. 

The comparison of the reaction of the adaptive loop to wind gusts for 
the various flight cases and aerodynamic conditions is shown in Table IX. 

A one degree change in the angle of attack corresponds to wind velocities 
greater than would be expected by the missile during the flight. The per 
cent change in for all flight cases and conditions is small, within 11 

per cent. The values of are returned to the undisturbed magnitudes 
while the gust is still being experienced for all flight cases and condi- 
tions in small inter»/als of time. The per cent change in and the 
time interval to return to conditions that existed before the disturbance 
increase as the airframe becomes mote flexible. 

The Reference signal is set so that if the dither signal source fails 
the value of is increased to insure that the system will not become 

unstable. The comparison of the rates of change of for the various 

flight cases and aerodynamic conditions when a dither failure occurs is 
shown in Table X. The rate of increase of decreases as the airframe 
becomes more flexible. 

The results of the investigation contained in this Chapter show that 
for the computer simulation of the dither system the performance of the 
adaptive loop deteriorated when the airframe became flexible. The adaptive 
gain is varied by actuator signals and wind gusts while the ideal system 
would have no variation of gain when experiencing disturbances. The 
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■niiLE IX 



.iiparisoii the Reaction of the Adaptive Loop to Wind Gusts 



1 Flight Case 
Aerodynamic 
Condition 3 


Magnitude of 
Change in 
Angle of 
Attack in 
Degrees 


Change in 


Per Cent 
Change in 


Time (seconds) 
to Return to 
Original 


Gain 
Setting 
in Adan- 
tive Loop 


t = 20 


-4.0 


0.0 


0.0 


0 


.025 




+4.0 


C .0 


0.0 


0 


.025 


t = 75 


-1.0 


+0.04 


4.0 


6.4 


.025 




+ 1.0 


-0.1 


10.0 


6.4 


.025 


t = 101.9 


-1.0 


-0.4 


1.7 


2.4 


0.2 




+ 1.0 


-0.4 


1.7 


2.4 


0.2 


t = 200 


-3.0 


-0.6 


10.7 


11.0 


0.1 




+ 3.0 


-0.6 


10.7 


14.6 


0.1 


< Flight Case 












Aerodynamic 












Condition 1 












t = 20 


-1.0 


0.0 


0.0 


0 


.025 




+ 1.0 


0.0 


0.0 


0 


.025 


t = 75 


-1.0 


+0.02 


1.8 ! 


4.0 


.025 




+ 1.0 


-0.09 


8.2 


4.4 


.025 


t = 101.9 


-1.0 


0.0 


0.0 


0 


0.2 




+ 1.0 


+0. 1 


4.5 


1.2 


0.2 


t = 200 


-3.0 


I 


0.0 


0 


0. 1 




+ 3.0 


1 0.0 


0.0 


0 


0.1 






* 1 


} 


( 
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lABLE X 



Comparison of Rates of Change of the Adaptive Gain for Dither Failures 



FligV»t Case 

Aerodynamic Condition 3 


Increase Per Second 


t = 20 


1.2 


t = 75 


1.12 


t = 101.9 


5.6 


t = 200 


1.8 


Flight Case 

Aerodynamic Condition 1 




t = 20 


1.2 


t = 75 


3.8 


t = 101.9 


6.2 


t = 200 


3.2 



i 






\ 



r 









£ 



I 



I 






variations of gain for wind gusts were small and were corrected during 



the gusts in a relatively short period of time. The variations of gain 
for actuator signals were large for two flight cases and small for the 
remaining flight cases. The variations of gain were corrected during 
the signal in relatively shor* periods of time. The rates of change of 
adaptive gain were satisfactoi ly based on the criterion of a constant 
rate of change of gain required between flight cases. By optimizing 
the gain of the self adaptive loop and properly shaping the adaptive 
gain pot the rates of change of gain could be increased, the variation 
in the adaptive gain due to actuator and wind gusts could be greatly 
reduced, and the time required to correct for disturbances could be 
reduced . 

3.5 Self Adaptive Loop . 

The importance of optimizing the gain of the adaptive loop was shown 
in Section 3.4. For small perturbations the loop is linearized in order 
that a qualitative investigation may be conducted (8). The linearized 
loop with the Reference value as the input is shown in Figure 3.1. The 
Reference value is compared with the rectified filtered dither signal, 
and an error signal is fonned. The error signal drives the integrating 
servo with the output shaft positic^ Indicated by (j) . The adaptive gain 



Kj versus servo shaft position. The gain, , is then multiplied by 



the Closed Loop Response Amplitude versus curve at the Reference 

value for the particular flight under consideration. The resulting signal, 



0, is then fed back through the filter. The characteristic equation of the 




which is the slope of the curve of 
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figure 5.1. Linearized Sell' Adaptive Loop, 



t 





1. The root locus is shown in Figure 3.1. 



As the loop gain increases, the response oscillates. It was found dur- 
ing the analog computer study that if the gain of the adaptive loop was 
increased beyond a certain value for a particular flight case that the 

response would oscillate and as the gain was reduced the response 

became over damped. The observations of the study agree with the root 
locus . 

In order to maintain the adaptive loop gain constant at the optimum 
value that is determined, the product of the Adaptive Gain Pot transfer 
function and the Vehicle Dynamic transfer function must be a constant. If 



of that the system selects for the particular flight case, the value 

of Kv will be a constant equal to 1. The shape of the gain curve for the 
pot for aerodynamic condition 3 with a dither frequency of 30 radians per 
second is shown in Figure 3.2. 

In an actual design a fixed loop gain could be established by varying 
and Kp until a desirable time response was obtained for the servo out- 
put shaft, and the gain of the pot could be shaped as explained above. 

If the performance of the loop was not satisfactory as to the rate of gain 
change, the change in gain caused by disturbances, or the time required to 
correct the gain after a disturbance, compromises would have to be made. 

The slope of the gain curve could be varied to improve the performance. 

In varying the slope, however, Kv would vary and the adaptive loop gain would 
no longer be a constant for a flight. A varying Kv would vary the time 
response of the servo shaft, and the performance of the loop would be af- 
fected. A different dither frequency with its corresponding Reference 
value could vary the shape of the gain curve resulting in the satisfactory 
performance of the loop. Optimum setting of the loop gains and shape of the 



the pot 




is the reciprocal of 




at the value 
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for the pot could be determined by a trial and error procedure. 
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CfiAFTER 4 



CLOSED LOOP RESPONSE AMPLITUDE AS A CRITERION FOR CONSTANT DYNAMIC RESPONSE 
4. 1 Introduct ion . 

The object of this Chapter is to determine for the dither system the 
accuracy of evaluating the dynamic performance of a variable parameter con- 
trol system - a flight control system - by the amplitude of the closed loop 
response, A dither signal is used as an input to the system and is varied 
by the dynamic characteristics of the system. The dynamic characteristics 
are then evaluated by the comparison of the amplitude of the dither signal 
with a fixed level Reference signal which represents the desired dynamic 
characteristics. The optimum dynamic performance of the system under in- 
vestigation was defined in Section 2,3 as a damping ratio of the short 
period mode of 0,7 with a minimum variation of the natural frequency. 
Therefore, the accuracy of evaluating the dynamic performance by the ampli- 
tude of the closed loop response is based primarily on the accuracy of 
evaluating the damping ratio of 0.7. 

The Reference signal values for dither frequencies of 20 and 30 
radians per second were not determined for an average damping ratio of 0.7 
for the flight for the computer similetion discussed in Chapter 3. The 
Reference values used, however, resulted in approximately equal average 
damping ratios for the flights, but the maximum variation of the damping 
ratios was twice as large for the dither frequency of 30 as for the frequency 
of 20. The maximum variation of the damping ratios of a flight also differed 
with the flexibility of the airframe and with the consideration of variable 
fuel load. The effect of varying the dither frequency, of increasing the 
flexibility of the airframe, and of varying the center of gravity on the 
accuracy of the performance evaluation will be determined in this Chapter. 

The results of the investigation will permit the selection of the Reference 
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signal value and the dither frequency^ or frequencies^ which result in 
the optimum variation of the damping ratio and natural frequency of the 
short per iod mode from the desired system dynamic performance a 

The investigation was conducted principally from an empirical view- 
point although some analytical verification is included in Appendix C« 
Aerodynamic Conditions 1, 2, 3, and 4 and Flight Cases t = 20, 75, and 
101«9, which are described in Section 2,4, are included in the investigation. 
Flight Case t = 200 was not included because the sensors are located in a 
different location than the sensors for the other cases. The performance of 
flights evaluated from cases that are measuring the reactions of the systen 
at different positions in the airframe was not considered a true indication 
of the capabilities of the dither system. Therefore, the performance of 
flights are evaluated in this investigation with only flight cases that 
are measuring reactions at the same airframe position. Application of 
the results to a general variable parameter system is discussed in Section 
4,4, 

4,2 Procedure , 

The investigation of the effect of varying the dither frequency on the 
accuracy of the closed loop response amplitude in determining the damping 
ratio of 0.7 was conducted using Cises t = 20, 75, and 101.9 as representa- 
tive of a typical flight. The dither frequency was varied in the range 
from 0 to 40 radians per second. The upper limit of 40 radians was select- 
ed because available actuators will not respond to higher frequencies. For 
various frequencies in the dither range of each of the four aerodynamic 
conditions Reference signal values were determined which resulted in an 
average damping ratio during the flight of 0,7. At the Reference value for 
a particular aerodynamic condition and dither frequency the maximum varia- 
tion of the damping ratios during the flight was determined. Plots of the 



54 



I 

I 

J 

I 

I 

I 

I 

I 

( 

I 

I 

I 

I 

I 



i 












I 



I 






i 

i 



I 

\ 




maximum damping ratio vat iat ion, , versus the dither frequency^ 1// , 

were then made for the four aerodynamic conditions. The procedure employ- 
ed in determining the Reference value and the maximum damping ratio varia- 
tion will be described by the use of a typical example, A dither frequency 
of 22 radians per second and aerodynamic condition 3 will be used for the 
example . 

Root loci for the three flight cases were determined for aerodynamic 
condition 3 using a Recomp II digital computer. The damping ratios, (f , 
of the short period modes associated with the adaptive gains, , were 

determined from the root loci, and plots of versus were made for the 

flight cases as shown in Figure 4,1. The amplitudes of the condition 3 

closed loop transfer functions were evaluated at the dither frequency of 
22 radians per second for various values of adaptive gain. The plots of 
closed loop response amplitude versus adaptive gain for the flight cases 
are shox^m in Figure 4.2, A graph of the damping ratio versus flight case 
was constructed as shown in Figure 4.3. The procedure used in establishing 

the closed loop response amplitude which would result in an average damping 

ratio during the flight of 0.7 was as follows: 

1. An amplitude for the closed loop response was selected, and 

the values of Kg corresponding to the amplitude were determined for the 
three flight cases from the Amplitude versus plots. 

2. The values of the damping ratios corresponding to the values 

of , which resulted from the selected amplitude, were determined from 

the cT versus plots. 

3. The values of damping ratios were plotted on the damping ratio 
versus flight case graph. 

4. If the average of the three damping ratios corresponding to 
the same closed loop response am.plitude was not 0,7, a different response 

amplitude was selected. Steps 1 through 4 were repeated until the average 
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damping ratio was 0.7o 

When the closed loop response amplitude which resulted in an average 
damping ratio of 0*7 was determined, the difference between the maximum 
and minimum values of the damping ratios for the three flight cases was 
taken as the maximum variation of the damping ratio, 9 during the 

flight for a dither frequency of 22 radians per second* The maximum 
variation and dither frequency were plotted on the versus 2^ graph, 

shown in Figure 4.4, for aerodynamic condition 3. The procedure was 
repeated for various frequencies for the various aerodynamic conditions 
to construct the versus l^plots for the dither frequency range. 

The optimum dynamic performance was defined as a damping ratio of 0.7 
with a minimum variation of the natural frequency. The procedure that has 
been described determined the variation of the damping ratios with varying 
dither frequencies while maintaining an average ratio of 0.7 during a flight. 
The procedure that will be described determined the variation of the natural 
frequencies of the short period modes with varying dither frequencies while 
maintaining the same average damping ratio as that used in determining the 
damping ratio variation. The reference for determining the damping ratio 
variation was the defined optimum performance ratio of 0.7. The reference 
for determining the natural frequency variation was selected as the average 
natural frequency of the short period modes of the three flight cases at the 
damping ratio of 0.7. The procedure employed in determining the natural 
frequency variation will be described by the use of the same example that 
was used for the damping ratio variation: namely, a dither frequency of 22 

radians per second and aerc dynamic condition 3. 

The natural frequencies of the short period modes at a damping ratio 
of 0.7 were determined from the root loci for the three flight cases. 
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g Figure 4.4. Daicping Fatio Variation vs Dither Frequency Ff- 
For Aerodynaiiiic Condition 1, 2, 'ana-4. : 
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The three frequencies were averaged to deteririine the reference natural 
frequency for condition 3. The associated adaptive gains^ ^ and 

natural frequencies, of the short period modes were determined from 

the root loci, and plots of Ijjn versus were made for the flight cases 

I as are shox^-n in Figure 4,5. The procedure used to determine the frequency 
variation was as follows: 

i 

1 1, The closed loop response amplitude which resulted in an 

I 

( average damping ratio of 0.7 for the flight was determined from the damping 
ratio versus flight case graph for the dither frequency of 22. 

2. The adaptive gains for the flight cases which correspond to 

the amplitude which resulted in the average damping ratio of 0.7 x^ere deter- 
mined from the Amplitude versus plots. 

3. The adaptive gains determined from the Amplitude versus 

j plots were used to enter the versus plots to determine the natural 
j frequencies of the three flight cases. 

' 4. The three natural frequencies were averaged. The reference 

I 

I natural frequency x^as then subtracted from the averaged frequency to 

determine the natural frequency variation for the dither frequency of 22 
radians per second. 

The natural frequency variation, A 7/A.^ , was plotted on the 

versus 7Z/ graph which is shown in Figure 4-6 for aerodynamic condition 3. 

The procedure was repeated for various frequencies of the various aerodynamic 

Ji 

conditions to construct the AZ^versus 7^ plot s for the dither frequency range. 
1 The plots of hi versus ~UJ and versus for the four aerodynamic 

, conditions will be used in Section 4.3 to determine the accuracy of the 

closed loop response amplitude as a criterion for constant dynamic response 
and the effect on the accuracy of varying the dither frequencies, increasing 
« the flexibility of the airframe, and varying the center of gravity. 
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^ ^ Results of the Linpiiical Invest Igat lon e 



One of the basic operations, which are discuss.d in Section 2^1^ 
of an adaptive loop is the continuous measaremerit of the system dynamic 
performance. The measuring method must have a regliglble effect on the 
system response signal. In the dither self adaptive flight control system 
the measuring mefhod involves a dither sig’nia^^ and the important system 
response signal is the short period mode responsoo The possible frequencies 
of the short period mode response are never greater than 1<,5 radians per 

• second for the cases and aerodynamic conditions considered. To insure 

'I 

" that the dither signal does not interfere with the short period mode re- 

j sponse the lower limit for the dither frequency range was set at 10 radians 

per second. The upper limit was previously set at 40 radians per second 
^ to conform to the characteristics of available actuators. 

I The effect of varying the dither frequency, increasing the flexibility 

of the airframe, and varying the center of gravity on the accuracy of evalu- 
ating the specified damping ratio is shown in Figure 4.4, The AS versus 
HJ plots are based on the requirement that an average damping ratio of 0.7 
I is maintained for the flight. The accuracy with which the damping ratio is 

j evaluated for aerodynamic conditions 1 and 2 is independent of the dither 

I' 

I, frequency in the range from 10 to 40 radians per second. The damping ratio 

|| for aerodynamic condition 1, which is the rigid airframe condition^ is very 

i 

I accurately evaluated with a variation of 0.06 between the high and low values 

I of the damping ratios during the flight. As the airframe becomes flexible 

\ by including the first bending mode, aerodynamic condition 2, the accuracy 

I decreases. The variation of the damping ratios for condition 2 is approxi- 

I 

I mately 0.15. As the flexibility of the airframe increases by including bend- 

ing modes 1, 2, and 3 but eliminating mode 2 by a bending mode cancellation 
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method, aerodynamic condition 3, the accuracy of evaluation the damping 
ratio depends on the dither frequency, A difference in the variation of 
the damping ratios for dither frequencies of 20 and 30 radians per second 
was determined in the computer simulation,, discussed in Chapter 3^ for 
this aerodynamic condition. The average variation of the damping ratios 
for the dither frequency range is approximately 0»24 with a minimum varia- 
tion at a dither frequency of 23 radians per second of approximately the 
same magnitude as that for a rigid airframe: namely j, 0o06. If bending 

modes 1, 2, and 3 with bending mode 2 eliminated by a frequency tracking 
and notch filter method and a varying center of gravity in the form of a 
slosh mode are included in the airframe equations^ the accuracy of evalua- 
ting the damping ratio is independent of the dither frequency from 12 to 
35 radians per second. The variation of the damping ratios for aerodynamic 
condition 4 in the 12 to 35 frequency range is approximately 0^435. For 
dither frequencies greater than 35 the variation decreases and is 0.375 at 
a frequency of 40 radians per second. The different methods used in eliminat- 
ing the second bending mode in aerodynamic conditions 3 and 4 resulted in 
large differences in the values of the real parts of the zeros associated 
with the third bending mode of the airframe transfer function as is discussed 
in Section 3.2. The zero for the third bending mode for condition 3 is close 
to the imaginary axis while the zero for condition 4 is not. 

In summary, the frequency in the dither range has no appreciable affect 
on the accuracy of evaluating the damping ratio for aerodynamic conditions 
1, 2, and 4. The accuracy depends greatly on the choice of the dither fre- 
quency for condition 3 with the magnitude of the variation varying from the 
small values obtained in condition 1 to the large values obtained in condi- 
tion 4. The accuracy decreases as the airframe becomes flexible and as a 
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varying center gravity Is Included. The large variations of condi- 
tion 4 could be reduced by using an average damping ratio of greater 
than 0.7 for the flight because for large damping ratios the amount of 
change of gain per unit change in damping ratio is greater than for small 
damping ratios as is shown in Figure 4.1. 

The effect of varying the dither frequency^ increasing the flexibility 
of the airframe, and varying the center of gravity on the natural fre* 
quency variation is shown in Figure 4.6. The variation of the natural 
frequencies for aerodynamic conditions 1 and 2 is independent of the dither 
frequency in the range from 10 to 40 radians per second. The variation for 
condition 1 is approximately 0.034 radians per second and for condition 2 
is approximately 0,123. The reference natural frequency for condition 1 
was 1.355 radians per second and for condition 2 was 1.185, The variation 
of the natural frequency depends on the dither frequency for aerodynamic 
conditions 3 and 4. The average variation of the natural frequency for 
condition 3 for dither frequencies in the dither range is approximately 
0.068 radians per second. At dither frequencies of approximately 23 and 
24 radians per second there are no variations. The average variation of 
the natural frequency for condition 4 is approximately 0,06 radians per 
second with a minimum variation of 0.028 at a dither frequency of 40 radians 
per second. The reference natural frequency for condition 3 was 1.225 
radians per second and for condition 4 was 1,51, 

In summary, the variation of the natural frequency does not depend on 
the dither frequency for aerodynamic conditions land 2 but does depend on 
the dither frequency for conditions 3 and 4. The variations were small for 
all conditions with 10 per cent of the corresponding reference value being 
the maximum variation determined. The effect of the flexibility of the 
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airfraiT’e and of a varying center of gravity on the variations of the 
natural frequency is not clearly indicated* The magnitudes of the average 
variations of the natural frequencies for the aerodynamic conditions varied 
from large to small in the following order: condition 2^ 3, 4, and 1. 

The investigation of the variations of the damping ratios and natural 
frequencies was conducted with the amplitude of the closed loop response 
selected to result in an average damping ratio during the flight of 0.7. 

Since the defined optimum performance of the system was a damping ratio of 
0.7 with a minimum variation of the natural frequency, the optimum dither 
frequency, or frequencies, for each of the aerodynamic conditions was 
determined on the basis of minimizing the variation in damping ratio and 
natural frequency. The error criterion used to determine the minimization 
of the variation in damping ratio and natural frequency was 
Plots of A5 versus are shown in Figures 4.7 through 4.10 for four 

aerodynamic conditions. The optimum dither frequencies selected for the 
aerodynamic conditions are shown in Table XI with the corresponding error 
criterion value, the variation in the damping ratio, the variation in the 
natural frequency, and the Reference signal value. Two dither frequencies 
with different error criterion values are listed for condition 4. The dither 
frequency of 40 has the lower criterion value, but the Amplitude versus 
curve for that frequency has a very small slope at the Reference value which 
would require a very accurate setting of the Reference signal. The dither 
frequency of 36 has a larger slope and would be the preferred frequency if 
the system was actually mechanized. 

The Reference signal values listed in Table XI are equal to the closed 
loop response amplitudes required for an average damping ratio of 0.7 for 
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TA.BIF Xi 



Opt inum Dither Frequencies and R. Values 



Aerodynamic 

Condition 


Opt imum 
Dither 
Frequency 
or 

Frequencies 
(Rad/SeCo ) 


Error 

Criterion 

Value 


Variation 

iv 

D-n-pi.-'e; 

Ratio 


Variat ion 
in 

Natural 
Frequency 
{Rad/Seco ) 


Reference 

Signal 

Value 


Condition 1 


10 








0.56 




26 


0.068 


0,06 


0,032 


0.485 




40 








0,42 


Condition 2 


10 








0,835 




12 


0,19 


C, 15 


0,118 


0,79 




26 








0,7 


Condition 3 


22 


0,105 


0,095 


0,047 


0,43 


Condition 4 


40 


0,375 


0,37 5 


0,027 


0,993 




36 


0,405 


0,405 


0,028 


0,975 
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tho The Is c xrC'-t if th*. - r y— ■. ^lue of 

the dither inpa^ si.'en^l r^d the g-iins :>f th- t . 1 ‘ -nd rectifier in the 

adaptive loop are equal to unity. If the ^alcis unity^ the 

Reference value listed must be corrected by being multiplied by the pro- 
duct of the half cycle average of the filter^, and the gain of the rectifier, 

4 , 4 Application of the Results to a General Case ^ 

The object of this section is to formulate a procedure from the empiri- 
cal results discussed in Section 4.3 and from the analytical verification 
contained in Appendix C to facilitate the rapid selection of the optimum 

dither frequency. The selection will be based a minimising the damping 

ratio variation. Once the dither frequency or area of frequencies is 
selected, the procedures of Section 4o2 can used for the frequency to 
determine the Reference value and the magnitude of the damping ratio varia- 
tion for the average damping ratio desired for the flight. 

The dither range is established by selecting the lo^er frequency limit 
five or more multiples greater than the maximum short period frequency to 
eliminate coupling. The maximum frequency of the short period mode can be 
approximated by the square root of the square of the real and imaginary 
parts of the zeros of the Inverse Models The upper limit of the range is 

fixed by the actuator requirements at 40 radians per second. 

The assumption is made that the control system will be investigated by 
a point study of the flight trajectory. From the root loci shown in Figure 
4.11, it Is noted that the open loop poles associated with the bending modes 
remain relatively stationary during a flight but that the ^pen loop zeros 
associated with bending modes move during the flight from smaller to larger 
values. If slosh modes are temporarily Ignored^ ch- effect that the varying 
values of the bending zeros have on the shore period mode is tc vary the 
shape of the short period characteristic. The bending zer^s must be in the 
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ii ^ t ■ • J ^ r.ppr'- effect 

on the shep-- ihe fi i;-'bt isc- wirb tb h the smallest 

imaginary part, which is the earliest flight Vw • ^asid^red after launch, 
results in a short period character 1st i : whi.h h-.s a larger adapative gain at 
a particular damping ratio than any of the other cases ^ The open loop poles 
and zeros associated with the short per iod appr ^ on ^ilsr vary during 

the flight, but it Is assumed that the vari^at i ''^os ar^ regligible compared 
to the variations of the bending zeros « The values of adaptive gain at a 
particular damping ratio decreases as the imaginary part of the bending 
zero increases. The adaptive gains for a particular damping ratio are, 
therefore, inversely proportional to the time in flight « The dither system 
maintains the amplitude of the closed loop response C'^nstant daring a flight. 
A constant amplitude results when the product of the adaptive gain and the 
dynamic gain remains constant. The dynamic gain is the gain associated with 
the vehicle dynamic characteristics. In order to have a minimum variation in 
the damping ratios, the product of the dynamic gain and the adaptive gain 
corresponding to the desired damping ratio must be approximately constant 
during the flight. Since the adaptive gain varies inversely with the time 
in flight, the dynamic gain must vary directly with the time in flight for 
the product to be a constant, Tf the dither fre(quency is selected close 
to a bending zero of the initial flight casc.^ the dynamic gains will vary 
directly with the time in flight. The herding zero must have a value with 
a small real part and an imaginary part in the dither range. 

Tf the value of the bending zero has an imaginary part in the dither 
range but a large real part, the selection of the dither frequency for 
optimum damping ratio variation is independent of the position of the zero. 
The value of the dither frequency does not appreciably the dynamic 
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r ^ ' Mir.af.tlv V. » sr. • r d nipi g . l. . ■ -icdc During 

the fligh- ibe bending s do result in v,. ti. *ns t : ae ^idaptlve gains 

for the same damping ratio* Since the vai.. t ..r dicker frequency can 
not result in dynamic gains varying inversely t the short period gains^ 
large variations in the damping ratio may result o The requirement of deter- 
mining when the real part of the value of the bending zero is large or small 
can be eliminated by always selecting the of the dither frequency 

close to the imaginary part of the bending zerv: of the initial flight case. 

If bending zeros do not exist in the dither rangey the optimum damp- 
ing ratio variation is independent of the selection of the dither frequency. 
The value of the dither frequency can not vary the dynamic gain to result 
in adaptive gains which correspond to the values required for constant short 
period damping ratios. The variations In the adaptive gains of the short 
period modes for the same damping ratio increase as the number of bending 
modes increase; and therefore^ the magnitude of the damping ratio variation 
would be expected to increase as the number of bending modes increase. 

Slosh modes may have an important effect on the dynamic performance 
of the system for some airframes^ Slosh modes would have a negligible 
effect during the initial and final phases cf the flight because the fuel 
tanks are then approximately full or empty. The slosh modes would affect 
the middle phase of the flight. If the slosh mode frequencies are close 
to the short period frequencies^ the shape of the short period characteristic 
could be appreciably varied for the middle flight caseSo The shape of the 
characteristic would not be affected for the initial and final flight cases. 
The largest adaptive gain for a particular damping ratio would be for the 
middle flight case which is most affected by the slosh modes. The dynamic 
gain would have to be reduced by the selection of a dither frequency for 
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• ' >'/C'.vi-ed. A dither 



frequency at j valuie (''^^.',11 t th.. ^ i . L'^'C'Crn^ of the 

middle case will rech.Ke. the dy-vcimic g'^in if the real part of the zero is 
small. The dither frequency selected^ however should be between the values 
of the imaginary parts of the bending aero of th^- initial case and the bend- 
ing zero of the middle case to result a small damping ratio variation 

throughout the entire flight® The exact frequency in the area would have to 
be determined by the procedure in Section 4® 2® 

In summary, the selection of the optimum dither frequency is as follows! 

1. If no bending zeros exist in the dither range^ the value of the 
dither frequency is not restricted, 

2. If there are no important slosh modes in the system but if 

there are bending zeros in the dither rangey the optimum frequency is close 
to the imaginary part of a zero associated with the initial flight case, 

3. If there are important slosh modes in the system and there are 

bending zeros in the dither range ^ the optimum frequency area is between the 
value of the imaginary part of an initial flight case zero and the imaginary 
part of the zero of the middle case v;hich is most affected by the slosh modes. 

The dither frequency must not assume a value that is close to the value 
of the imaginary part of a pole which is adjacent to the imaginary axis be** 
cause the adaptive gain change per unit amplitude of the closed loop re- 
sponse is large, A large gain to amplitude ratio requires a very accurate 
setting of the Reference value. 
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CHAPiER 



CONCLUSIONS 

The conclusions from the results of the investigation are as follows: 

1, The dither adaptive loop employs a single parameter to control 
basically a single dynamic characteristic of the system^ The parameter of 
the loop is the closed loop response amplitude and the. characteristic con- 
trolled is the damping ratio. 

2, The components of the loop are simple to mechanizes, but a prior 
knowledge of the airframe dynamics are required to select the dither fre- 
quency, Reference value, and the shape of the gain curve for the adaptive 
gain pot. 

3, The adaptive loop will adjust the adaptive gain to partially 
correct for the varying dynamic characteristics of a flexible airframe if 
the dither frequency has no coupling with a bending mode^ The accuracy 

of the correction in maintaining a constant damping ratio is the inherent 
accuracy of the Amplitude Closed Loop Response Criterion. 

4, The adaptive loop makes satisfactory corrections of the adap- 
tive gain for disturbances. A further investigation into the gain curve 
shape for the adaptive gain pot and into the optimum lo:p gain is required 
to verify that an actuator signal will not cause an unsatisfactory variation 
in the adaptive gain. 

5, The Amplitude Closed Loop Response Criterion does not result 
in a constant damping ratio. The variation of the damping ratios increases 
as the desired average damping ratio decreases^) as the airframe flexibility 
Increases, and as the variation of the center of gravity increases. If 
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£requericics in the cMthrr range greatly affect the amplitt.de of the vehicle 
dynamic response, the variations of the damping ratios can be appreciably 
reduced by the proper selection of the dither frequencyo 

6* When the damping ratio variation varies with the dither fre- 
quency, the minimum natural frequency variation occurs at approximately the 
same dither frequency as the minimum damping ratio variation* 
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The leiodyTicirnir x^^quations for the missile were cbtained from the 
Aetonetics Division of North American Aviation. The equations are two 
degree of freedom shrt period node approximations with three bending 
modes and one fuel slosh, mode Included, 
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-.00704 


-.0052 


cr 


0 


-.0482 


-.0125 


-.0013 




0 


-..18 


-.43 


10.24 


Er 


0 


.0469 


.0.24 


.015 


Er 


0 


1.27 


.584 


12.7 


Gr 


-.832 


-.906 


-.852 


-1.49 


V 


-1730 


-1734 


-1732 


-5351.8 


V 


-.0182 


.476 


.18 


0 


“V 


.0374 


.085 


.181 


3,52 


V 


.0288 


.0655 


-.14 


.018 


Lr 


-1.57 


-1.78 


-1.79 


-2.18 


A; 


-2.3 


-2.3 


-2.3 


0 


Xsl 


3.4 


3.4 


3.4 


3.4 




-6.95 


-6.95 


-6.95 


0 



The magnitudes of the remaining coefficients were not available but are 
included in the transfer functions. 
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Infomat ion 



' , f t t -.M ions 0 and D< will be shown. 

4 f » » • • t t 

L( ,1 't ’in .- tRi j f nation of the d, ^ , dj , Ai , AloX 

equal 1-Jiv- is contained in i\e Terence 10« 



and Atrt 



Fron ^^ex^ton*.s second law of motion. 
^P/ - J ^ 



^ 13:^ 



(2) 



+-^^(S7xA ) 

i.// r 



where ^ Fy , ^ F^ are the forces in the x, y, z directions of a 

right hand system of Cartesian axes fixed in space; where ^L, ^ M, 
are the nenrents about the x, y, and z axes; where ^ ^rid are 

the moments of momentum about x^ ^ which are fixed to 

the airframe; \^here J^(w X X h) , JL{U/ )t^]r\ ^ arise from 

the angular velocity 'HJ of the Xj , , and zj axes with respect to 

the X, y, z axes; v/here , ^/y , accelerations in the 

X, y, and z directions; and where m is the mass of the airframe. 

The absolute acceleration of a body with respect to x^ z axes is 



where V Is the instantaneous linear velocity and 24/ is the angular velocity, 

() + ✓+• vi/ 

^/✓ = (Qw/-Rv)^ r(ftu-Pw)^t(PV-q>u)J' 

where F, Q, and R are the angular velocities along the Xj y, and z axes 
respectively and U, V, and VI are the linear velocities along the x^ and 
z axes. Thi components of the acceleration are then 

^ U (pW^ R X 
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\ 



( 

I 

1 



r 



' -- V + RU-Pv/ 
o.y^W rPV-QU 

Substituting the comp<;nents of acceleration into Equations 1 and consider- 
ing the mass as a constant the equations become 

(utQW~RV) 

£Fy ~ ^yu(v -r RU - Pw) 

£~Fz ~ A>^-Cw P \/ " c^u ^ 



The components of the moment of momentum are calculated from summing the 
moments of the velocity vectors about each axis and multiplying by the mass ci/w, 

hy " CpSyy-- PXyy^ ""PX ^ 

where r and T-jOy - 

XY and XZ are assumed planes of symmetry, and therefore 



(^) 



l^y 






7^ 



( 5 ) 



(6) 



- 0. Equations 4 become; 

h<^ “Pi 

h-y ~^1yy 

hy' 





- I-'y/y 

H X ^ * 

ilixh --(h^.Q-hyR)^+(h^R-k^P)/+{hyP-f,^Q)X 

Equations 2 with Equations 5 and 6 substituted become 



Since 



86 



( 7 ) 



" L - 

Ql.yAy t "'^'^'y^ 

2 . a / - i\~l^,/^jr P(p(Tyy — • 

The forces indicated by the left terms of Equations 3 represent the 
summation of aerodynamic, thrust, and gravity forces. The gravity forces 
are shown below and their derivation can be found in Reference 9. 



(8) 



~ C" 



Oo)cr^0 



Cck2./^ 



.) 






^ ^o) OfP^ ^ 

* C*** o ) Co^ (p^"t 

^ Cc^Sf,<iAri^ 

Ca%i^Sa C<r<i^^(>)(c^ & 

P^/y'Z. /yry'^.-oA^y 0 CtrOj (f (j^) -h 

{ Cr^y' 00 (^0 0 (^ ~ Cr^ ^ ■<)-C'r\y(f^ -f- 



( 



y7h^ 






CyrK^^ 



0(y Csr^ cp^ • 



where m is the mass cf the body; g is the acceleration of gravity; e, <pj 
and y are the Eulerian Angles for a ZYX system and Qq^ and y^j are 



the steady flight values, 

/ 



Defining ^ ^ F y ^ and ^ as the sum of the 

aerodynamic and thrust forces. Equations 3 can be written as 



( 9 ) 



= 

£Fy = 

£pL- 

A 



(u-t Qw^ rv^ 

^n^(v 'hpi/-' Pw) ~~F^ 
/m--(VV^ ~hPV ~~ (pU } f 
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F.acii of thi; 1 ■ni;aneous velocity components can be written as the 



sum of a steady state value and a change in velocity. 



U ^ Uo f 

V ^ Vo f /V^ 

VV = V/oYuT 

P ^ Po 

Cp r QoY^ 

P ~ Pot 

^ 

c£^ 



In order to linearize the force and moment equations it is assumed 

that 

1) the velocity changes are small and that the products and squares 
of the velocity changes are negligible compared to the velocity changes 
themselves 

2) the disturbance angles are small and the sine of the angles 
can be set equal to the angles and the cosine of the angles set equal to 
one. The product of the angles is zero. 

Equations 7 and 9 can now be written as: 



(10) 



(H) 



-t QoVVo Qo/t^~~PoVo - /?o /tr'- i' 

~ ■/■ Oofto f" Uo-Ztz "t RcMz" R) WC " 

~ ~t foVo+ Pa AT' + " (^oUo "■ (pOyOO^ '•'Uo^ f" 

p' Bo) B Bo 

£i + ( Qo/?0 ■/'^po 
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^ M ^ j y y ^ G n f* t 

X ^ ~ I''^ ' ^ ^ V *" t ^'^o) L^^iAj '“ "T J ^ 

^ ^ 

The equations which relate angular velocities about the x, y, z axes 
to the Eulerian Angles 0y </>. ^ are 

Pc 

• « 

Cp - 0 Cd^a- iP xJtW p Ca^ 0 

R ' ^ OerO' p Ca<U 0 — 0 (p , 

When linearized the equations become 



(12) 



P-- {/; 

9= © 
R: 9 



The steady flight condition is defined as all velocity components equal 
to zero except and « Using this definition and Equations 12, Equa- 
tions 10 and 11 become 



(13) 



- />>1~ -h Wo^'' ir ^{c<rO &o) 0^ 

JLaT •+ Uo A^ ““VV<9^ ^ ^ e.)fy(^e,)cf>2 

£ = ^Lui^^Uo ^ t orfOo^ 



(14) 

£l = CpX^^ 

£ M C Qjyy 

£t/- 

The aerodynamic forces and moments can be represented by a Taylor 
series expansion* Because of the planes of symmetry and the steady flight 
conditions ^ functions of only q, u, and w of the 

velocity components and ^F.v; , and N are functions of only p, v, and r. 



ihi. !S , 'I v^ice 9, Iht? other variables affecting the forces 

and I ■ . Ls at ‘ the angle of thrust deflection, the airframe, flexibility, and 
the cVianging mass of the body. Since the changes ate assumed small the 
second derivative and above for the velocity and bending terms are consider- 
ed negligible and the first derivative and above is considered negligible for 
the change in mass. The resulting equations are of the form 



v;here d^ ^ , and d^ are the first, second, and third bending modes; 



the changes due to the change in the mass of the fuels <, 

The thrust axis is assumed to pass through the center of gravity and 
hence does not contribute to the moment equations. The thrust force equa- 
tions are 



05 ) 




Sr is the angle of thrust deflection; and AI , AlOX ? and are 





is the angle between the X axis and thrust line. 



+AT' 



and 




Therefore 




( 16 ) 
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To determine the e^^dations of steady flight substitute Equations 16 
and 15 into Equations 13 and 14^ substitute the values of the steady flight 
condition, and set the velocity changes equal to zero. The resulting equa- 
tions are: 



(17) 



'fyo O 
Lo 



Mo - O 
h/o 'Z O ^ 



Wlien stability axes are used as the reference system and quasi-steady 
flow conditions are assumed^ all terms containing disappear and all 

aerodynamic partial derivatives with respect to rates of change of velocit- 
ies and with respect to rate of change of angle of thrust deflection are 
eliminated. For the Two Degree of Freedom Short Period Mode Approximation 
is set equal to zero and only the equations ^ and are 

involved. Substituting Equations 16 and 15 into Equations 13 and 14 and 
subtracting Equations 17 gives the following equations of motion for the 
disturbed body: 

Oc ^ f Q o)B = +• fdtJt 'f' 

-h F -f a4^ ■l' Fd^3 dis i' 



(18) ^ • 

^ -f t M^tSt ^ ^ 

Mctl t M/i /i AT -h MmqxALOX i- 
AIH 



where 




J- 




and M <y 

/J 
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ibc ciP: ' 



f 



1 



o< 



d< fined as the angle between the wing 



chord line a .u t hi r. wind. 

/v. /\0<- U - ~Uo /O^-O 

(19) Th. ref. 'c S' [)oO{- 

Equations 18v » :h quations 12 and 19 substituted become: 

u a ^ M 9 , j) a t (^) 4 ^ 
( ife’) c// f (^r) /3^2 -t C ^4f) (£a f- (^)Js t 






e =(m )o< -t (ai^) Ot(Mdj)<^i '/' -h(Mdi) 

-h )d^. t f (M Jj) ^3 1 Cm^x)/xT 

■^Oi LOA )ALoX H'^alh)ALH. 

The coefficients for any instant are constants. Equations 20 agree in 
form with the given equations for oL and © # 

The Transfer Functions which were obtained from the Aerodynamic Equa- 
tions are ar follows: 

1. Aerodynamic Condition 1. 



So 


Gp - 


/. ■f. 0 


7S 


Qp ^ 


A • 




1 II -ti oAo'i 9 




Bp - 


S.o^ (^i-.D3lH) 






- Ado 


Bp - 




<^T 


(^->co 7H)C<!'-7 AS)(^iJ,3A') 
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lyricij.ic 



/r~c3o 






&p _ 

ir 



n 2. 



y ^ 6 ^ 0,34^ ^ rj ) 

y^(y<^l 0 ) Co i~, o $> t C> ) 






^ - (^.2^C^-h.o(>rJ(^.o^^^^s: ! ) 

,oa3)U'-> 77 )(^ ■*'- ^9)U-i-> -^' r. «2 ) 



/'ph)> ? 



■5r 



9. J f . t>32)(-i>'t . qSJX. ^J>lC^.1 

o)) (a^ " . ^ T*"/ 79/^ j', 0 ’^1 B> I7^ 



^OO 



|L= ^ ?: 

'^7' (<i.'',oc7) 7 ) 



3, Aerodynamic Core it ion 3; 






7£ c 3 p - / 9 . / r<w- « 0 ^ 6 ^ Ci 2 - '^<6 67 1 ^ V. 98)C^-t > ! ^ iJJT 

~ , Oe^ 3^ Cc' "^1 7 7^ (^ y*, S’ }o3^^ 8> ^9^^'hi 



y - )ch ? z ^!> ioC-oi io2^)(’>i> i‘< i^oS’ -j!^A7, OS"^ t 97^ 

(a^-,c l)Co^. C>23jC>i.7H3^(o-r- a B7 ui.iUj 
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^ ^ 39)(^ +.3i9 jxr. 73 . 3 ) 

^ 7 - £^f^np^.<^8rr^'u6.7)4t,7ift^7s,ij 



. Aerodynamic Condition 4; 






( .0<3 ? - y lo^)(f> - /A Ih'^^.tj^ .^.l< k). 

i^i , (9 / /,9l>)i^-t<olS' t^' /, S'!) Cs^j- 



7^~7S - l^‘7U9'.ol)Ut.oe7^^-A90(^f.ofiAt^:7.j}p)/!^t,o^?tJ^.^O 

^ (^tAS'J U-.ol)U- 

( ^ 9, c 99 ^J^9. ^7^-^ •^io,A -ity* r3^, (? - ^ 

u,t\ oi> t !o3 1^ 2, 



^ 3p ~ JS , i^T-^ r>cH)C>'’9< 3t j^c^9,^9s£^2 (^^ ^ j 

r^r . ^34 jy. X? ^ ) la 1 ^. ^ ^ rj.^ .i^)C^t.o3s^ t^a> ^ ll. 

(^f .03 C t^' 9 A ^ H) (]ar f 7o8?f^. B. iS) U ^ . oJH t ^ 3. HSi) 



y^ ~ f^O'~> ^ ^ - f?. 2 S' {. .^9 , o 3^ ')C<!i^f^777^ei I J7V^ye^» 2j(d'i‘> 3 A 

cSr ^ (■<>■ -.CCU)(^^. ISStj^ J> 7^X(^-h.9l tJ,S<3 1)<^9, ±^c3. 77/ 

f . HS8- rx9^, 7^)U f-. 7rj^> U 73. j ) 

(^f-, ^y/ t^;?Z r tp HLl){or,79S 
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APPENDIX B 
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Figure B.l, Computer Simulation Diae^ram. 
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Figure B.l. Continued. 
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Figure B.l. Continued. 
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Figure B.l. Continued, 
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Jt’xgure Disturbaiice Heaction of system for i’lie,ht Case t=75 











FitjOire B.^(b). Continued 










Figure B.3(c). Disturbance Reaction of System for Flight Case t»101.9 
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FifcUre B. 5(d), Disturbance Reaction of System for i'li^dt Case t=200 




i'igure B.3(a). Continued 



ANAlYil-^AL Vinr^irArrOK OF i.tL ErFiHa.Ai 

The results of the variation of damping ratios >?ith dither frequency 
obtained by empir: ?A methods in Chapter 4 ^/il! be analyzed by analytical 
methods in this Appendix* 

For the analysis the root loci for the flipht cases of a particular 
aerodynamic cond tlon will be assumed to have the same shape and system 
gain values. The transfer functions of the vehicle dynamics for a parti- 
cular aerodynamic condition do not vary to the extent that the use of the 
assumptior to determine general trends would be preventedo A particular 
syste'^ gviin^ K, on the root loci for the three flight cases of a condition 
xvould determine the same damring ratio^^^o The system gain is the produc: 
of the dynamic gain, , and the adaptive gain^ K ^ o The dynamic gain 

Is the gain associated with the vehicle transfer function^ Using condition 
3 as an example, the dynar^ic galr 4 S for the flight cases are approximately 
equal to the following values: 



Case l 20 




16,05 




Case t ^ 75 


- 


19.1 




Case t 101,9 


- 


21,6 




Since a particular damping ratio has the same 


va 1 UP 


of K for the three 


cases, the values of for the cases 


are . 


K 




Jo 




TbTUT 






~ 


19.1 

K 




K^^ = /o/,9 




21,6 




and t - 20 > t ^ 75 > 


t - 


101,9 


for a particular ^ 



A plot of ^ versus : ^ Is of the follcwing form„ 
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K,$ 



The trends shown in tVie S versus plot are valid for any of the four 

aerodynamic conditions. 

The closed loop transfer function for the system is 

e _ 

COMM/i/VD Si&t^AL 

when vehicle dynamic transfer function. When the dither freq- 

uency is large compared to the poles and zeros of the vehicle dynamic trans- 



fer f mction, the transfer function is approximated by 



_a£ 



St 

This approximation is justified for the dither range^ 10 to 40 radians per 
second, for conditions 1 and 2. The closed loop transfer function becomes 



© 



CaAlMAf/D SIGML 






1 + 



Since the dither system maintains the amplitude of the closed loop response 
constant , 



Ks, N'U, 


- ks^MsA 




i Ms, 






and for Aerodynauuc uouaition 1 
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k 



k^, Isis' _ - ks^ S,o(> 

I -h ks, /. I i iUiSoi^ 



and for Aerodynamic Condition 2 

Ki..s:i^ - u.r - 9 . .3 

li-KiiS.iS \+Kij6.7S' H-«i3 9,A 

The relationships are independent of the dither frequency^^ and constant 
damping ratio variation would be expected in the dither range. This agrees 
with the empirical results for conditions 1 and 2 as are shown in Figure 
4.4o If is assigned the value 1 for both condition 1 and 2^ then for 

condition 1 



2 = 1*67 0«516 

and for condition 2 

1 - '-3 K^3 = 0-725 

For both conditions, 2 ^ 3 trend necessary 

for small damping ratiw variation as determined from the cT versus plot. 

For condition 3 the dither frequency can be close to poles and zeros 
of the vehicle dynamic transfer function. If the dither frequency is close 
to a pole of the transfer function^ then -jr- is very large and 



e 



and 

A 



COtAfy^Ar/O ~ + 

The amplitude of the closed loop response is independent of the adaptive 
gain; and therefore^ could assume a infinite number of values with 

an infinite number of damping ratio variations. 

If the dither frequency is close to a zero of Case t = 20^, then 



0; 






< 









.Or, 



St^z hi, 7 
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variation. If the 



:j anpljfuut of The closeu loop sp'>rtse t = 20 

t 101.^ which is the trend for a small damping ratio 
dither frequency is close to a zero of Case t = 101,9, 



then 

9 p \ Of \ &F , 

In order to have the same amplitude of the closed loop response Ks t = 20<; 

Kx r = 7^ = 101.9 which is the trend which would tiosuit in a 

large damping ratio variation. If the dither frequency is close to a zero 
of Case t = 75, the magnitude of the damping ratio variation would vary 
between the magnitudes when the frequency is close to a zero of Case t = 20 
and when the frequency is close to a zero of Case t = 101. 9<> 

The damping ratio variation curve for condition 3 should^ therefore, 
have an undetermined damping ratio variation at the pole at approximately 
40 radians, have a small value for the variation at the zero of case t 20 
at approximately 24 radians; and have a large value for the variation at 
the zero of case t = 101.9 at approximately 27 radians. This agrees with 
the condition 3 curve shown in Figure 4.4. 

In order to determine the dither frequency that the minimum variation 
in damping ratio occurs, the ideal expression for no variation which is 
of the form 




Jr 5f 



sr 




_J_ 3 =: L-, 

6 ! -j. I , a. 

' s n /I ~ j ^ ^ 



<Sr 






c^r 



i 






^'1' 

&p 



Sr 
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wht / ♦. 



n - I 

3 = KJ t 
C . t 






75 



> / 
go 
f- g o 



101.9 —‘^g— 

go 



and for a damping ratio of 0,7 
t - 20 - 0.86 

Kc t = 75 = 0.26 

C 7 

t 101.9 0.175 

Is evaluated for various dither frequencies close to the frequency of 
the zero for case c = 20, and the smallest amplitude variation between the 
three cases should indicate the optimum frequency. A dither frequency of 
23 radians per second resulted in the smallest variation in the amplitude. 

For aerodynamic condition h there are no zeros or poles of the vehicle 
dynamic transfer function which are close to the imaginary axis in the 
dither range. Varying the dither frequency does not result in wide varia- 
tions in the values of the vehicle transfer functions, and the closed loop 
transfer function can be approximated by 



0 

COMMAh/O SIGt/AL 



1+ KiMs 



A constant damping ratio variation would be expected and agrees with the 
curve for condition 4 in Figure 4.4. Since the values of the vehicle 
transfer functions can not be varied appreciably, the variation in the 
adaptive gain required for a constant damping ratio can not be realized 
and a large variation in the damping ratio results. 
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